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' Anisotropic systems of quarks and gluons, which at least for sufficiently short space-time intervals 

, can be treated as homogeneous and static, are considered. The gluon polarization tensor of such a 

'n— «' ' system is explicitly computed within the semiclassical kinetic and Hard Loop diagrammatic theories. 

The equivalence of the two approaches is demonstrated. The quark self energy is computed as well, 
j ■ and finally, the dispersion relations of quarks and gluons in the anisotropic medium are discussed. 
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I. INTRODUCTION 



> 

\^ • The state of equilibrium being static and homogeneous is sometimes anisotropic. This may happen when the system 
of quantum fields, which is of interest here, is under influence of an external force. A relativistic plasma, which is 
anisotropic due to a magnetic field, often occurs in astrophysical situations as the early Universe or Supernovae [0. 
Anisotropic states are also common for systems which are out of equilibrium. Sometimes such states can be treated 
as static and homogeneous, but only for sufficiently short time and space intervals. How short the intervals should be 
depends on the specific problem under consideration. 

The parton system generated at the early stage of ultrarelativistic heavy ion collisions at RHIC or LHC is of 
particular interest for us. The parton momentum distribution is not istotropic but strongly elongated along the beam 
O I Therefore, specific color fluctuations, instead of being damped, can exponentially grow and noticeably influence 

O the temporal evolution of the system. In a series of papers of one of us |^,^ it has been argued that there are indeed 
' very fast unstable plasma modes in such a parton system. The stability analysis has been performed within 
the semiclassical transport theory of quarks and gluons Wm . Since the theory has been proven till now to be fully 



consistent with the QCD dynamics only for quasiequilibrium systems |^,^ one wonders to what extend the results 
rN from jj,^ are reliable. Thus, a QCD diagrammatic analysis is desirable. 

1 Perturbative approaches within the real time field theory provides a natural framework to study weakly interacting 
■ ■ ■ ' quantum field systems in and out of equilibrium. However, the naive perturbative expansion, when applied to gauge 
fields, suffers from various singularities and some physical quantities are even gauge dependent. These problems have 
been partly resolved for equilibrium systems by using an effective perturbative expansion where the so-called Hard 
Thermal Loops are resummed |p^ . The Hard Thermal Loop resummation technique within the finite-temperature 
QCD has been shown to be equivalent to the approach based on the classical transport equations, where color 
is treated as a classical variable, or on the semiclassical j|] one, where the color degrees of freedom emerge due to 
the matrix structure of the parton distribution function. The Hard Thermal Loop approach has been generalized to 
nonequilibrium systems, but only very specific forms of deviations from the equilibrium have been discussed so far: 
systems out of chemical equilibrium, which are important in the context of heavy-ion collisions [ p^ , and such where 
the momentum distribution is isotropic but not of the Bose-Einstein or Fermi-Dirac form ||r^-[l5[| . As observed in , 
the Hard Thermal Loop approach can be applied to any momentum distribution of hard particles which is static and 
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homogeneous. This is evident when the Hard Thermal Loop effective action is derived within the transport theory 
. The term 'thermal' is then rather misleading and for this reason we shall omit it in the following. 
In this paper we discuss the applicability of the Hard Loop technique for systems with anisotropic momentum 
distributions. The technique has been earlier applied to the equilibrium QED plasma in a magnetic field |]l7| . Our 
aim is to consider a general situation with an arbitrary momentum distribution. We analyse the problem from the 
point of view of the transport theory and the diagrammatic approach. Using the scmiclassical kinetic equations we 
derive the Hard Loop induced current paying much attention to the gauge aspects of the procedure. We also explicitly 
demonstrate that the gluon polarization tensors found by means of the two approaches are identical. In this way, the 
applicability of the kinetic theory beyond the equilibrium is substantiated and more specifically, the reliability of the 
results from is shown. 

The Hard Loop diagrammatic technique has the advantage over the semiclassical transport theory approach that 
it can be naturally extended to fcrmionic self energies and to higher order diagrams beyond the semiclassical approx- 
imation. In this way the dispersion relations of quarks and other observables of the quark-gluon plasma, such as the 
energy loss of energetic partons, transport coefhcients, or photon and dilepton production rates can be calculated 
systematically in the case of anisotropic distributions. We take a first step in this direction computing the quark self 
energy for an arbitrary momentum distribution. 

The self energy controls the particle dispersion relation which provides an essential dynamical information about 
the system. We discuss therefore the general dispersion relation of gluons (plasmons) and quarks in the anisotropic 
quark-gluon plasma. Finally, we briefly consider possible applications of the formalism developed in this paper. 



II. TRANSPORT THEORY APPROACH 



In this section we first introduce the semiclassical transport theory of quarks and gluons ||6|,|7| . Then, applying the 
linear response method, the Hard Loop induced current is derived. Finally, we compute the gluon polarization. 



A. Transport equations 



The distribution function of hard (anti-)quarks Q{p,x) ((5(p,x)) is a hermitian Nc x Nc matrix in color space 
(for a SU{Nc) color group); x denotes the space-time quark coordinate and p its momentum. The four- momentum 
p = {E,p) is assumed to satisfy the mass-shell constraint. Since both quarks and gluons are treated as massless 
particles the constraint is = 0. We also mention here that the spin of quarks and gluons is taken into account as 
an internal degree of freedom. The distribution function transforms under local gauge transformation M as 

Q{p,x) ^ M{x)Q{p,x)M^x) . (1) 

The color indices are here and in the most cases below suppressed. The distribution function of hard gluons is a 
hermitian (A^^ ~ 1) x {^c ~ 1) matrix which transforms as 

G{p,x) ^ M{x)G{p,x)MHx) , (2) 

where 



Mab{x) = Tr[TaM{x)TbM\x)] 



with Ta, a = 1, iV^ - 1 being the SU{Nc) group generators in the fundamental representation. 
The color current is expressed in the fundamental representation as 



fix) = -g 



{2nf2E 



Q{p,x) - Q(p,x) ~ ■^Tt[Q{p,x)-Q{p,x)] 



2iTafabcGbcip,x) , 



(3) 



where g is the QCD coupling constant, fabc are the structure constants of the SU{Nc) group. 

The distribution functions of quarks and gluons are assumed to satisfy the following coUisionless transport equations: 
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d_l 
d_l 
d_l 



p''D^Q{p,x) + gp^' — -{F^Ax),Q{v,x)] = , 



p^'D^Qip.x) - gpt^ — -{F^,ix), Q(p,x)} = , 



p''V^,G{p,x) + gp'' — -{T^,Ax),G{p,x)} = , (4) 



where {..., ...} denotes the anicommutator; Z?^ and 2?^ are the covariant derivatives which act as 

Df, = df, - ig[Af,{x), ... ] , = 5,, - ig[Af,ix), ... ] , 

with and A^ being the mean-field or background four-potentials; 

A'^ (x) = A^, {x)Ta , ix) = -ifabcA^ (x) ; 

F^i, and J-f^,y are the mean-field stress tensors with a color index structure analogous to that of the four-potentials. 
The background field is generated by the color current (^) and the respective equation is 

D^F^^'^ix)^,f{x) . (5) 

We note that the set of transport equations (|[ ||) is covariant with respect to the gauge transformations (|l|, ||). 

B. Plasma color response 

We discuss here how the plasma, which is (on average) colorless, homogeneous and stationary, responds to small 
color fluctuations. The distribution functions are assumed to be of the form 

Qij (p, x) = n(p)(5y + SQij (p, x) , (6) 
Qij{p,x) = n{p)Sij + SQ,j{p,x) , 

Gab{P, X) = ng{p)Sab + SGab{P, x) , 

where the functions describing the deviation from the colorless state are assumed to be much smaller than the 
respective colorless functions. The same is assumed for the momentum gradients of these functions. The (anti-)quark 
and gluon distribution functions nljp), n(jp), rig(p), reduce in equilibrium to the Fermi-Dirac or Bose-Einstein form 
i.e. 

"^P^ ^ exp(|p|-M)/T+l ' 

-( ^ = ^^^^^— 
"^P^ exp(|pH-/i)/T+l ' 

"«^P^"exp(|p|/T)-l' 

where T and /i denote the temperature and quark chemical potential, respectively, while the factor of 2 occurs due 
to the spin degrees of freedom. The number of quark flavours is assumed to be equal to one. 
Substituting in (||) one gets 



j''ix) = -g J j^^^p''[sQ{p,x) - 6Q{p,x) - ^Tr[dQ{p,x) - 6Q{p,x)] (8) 

+ 2iTafabcSGbc{P,x) . 



As seen, the current occurs due to the deviation from the colorless state. Let us also observe here that not only 
(anti-)quarks but also gluons contribute to the current (^). Thus, the current is of essentially non-Abelian nature. 

Now, we substitute the distribution functions (||) into the transport equations (||). Assuming that the stress tensor 
is of the same order as SQ, SQ or 6G and linearizing the equations with respect to 6Q, SQ and SG we get 
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P^V^6G{p,x)^~gp^T^,ix 



dn{p) 

dpi, 
dnjp) 

dpi, 
drigip) 

dpu 



(9) 



We keep here the covariant derivatives to maintain the gauge covariance of the equations. 

To solve the equations such as Eqs. (^) one usually uses, see e.g. the gauge parallel transporter defined in the 
fundamental representation as 



U{x, y)='P exp -ig I dz^ 



where V denotes the ordering along the path from x to y. There is analogous formula of the gauge transporter U{x,y) 
in the adjoint representation. Using U and U one finds the solutions of Eqs. (0) as 



5Q{p, x)^-g j d^y Gp{x - y) U{x, y) p^F^^iy) U{y, x) 
SQip,x) = g f d^yGp{x-y) U{x,y) p''F^,{y) U{y,x) 



dn{p) 

dpu 
dn{p) 

dpu 

SG{p, x) = -5 J d^y Gp{x - y) U{x, y) p^^T^ ^{v^ ^) . 

where Gp{x) is the retarded Green's function which satisfies the equation 

p^d^'Gpix) = 6^^\x) 

and equals 

with t being the 0-th component of x (a;^ = (t,x)), and v denoting the parton velocity i.e. v = p/E. 

Substituting the solutions (|o|) in Eq. (||) one finds the color current of the gauge covariant form which reads 



(10) 



d^p 



{2Trf2E- 



p^'p^ / d^yGp{x-y) U{x,y) F^Av) U{y,x] 



dfjp) 
dpi, 



(11) 



where /(p) = n(p) + n{p) + 2Ncng{p). 

Now, we are going to perform the Fourier transform of the induced current (^l]). Before this step however, we 
neglect the terms which are not of the leading order in g. Then, the transporters U are approximated by unity and 
the stress tensor F^^, by di^iA^ — dyA^. Within such an approximation, the Fourier transformed induced current jll]), 
which is no longer gauge covariant, equals 



{2ny2E ^ dpx 



p'^ka + ^0^ 



A,{k) . 



(12) 



The induced current j^{k) can be expressed as 

fa{k)=UZ{k)Al{k), 

with n^"^ being the giuon polarization tensor. Transforming Eq. ( p^ to the adjoint representation one finds 



n-(.)^.7^ 



{2'Kf2E ^ dpx 



k^p 



p'^ka- + iO+ 



(13) 



It should be noted here that the polarization tensor is proportional to a unit matrix in the color space. 
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Now we are going to show that the polarization tensor is transversal i.e. A:^n^'^(A:) = 0. Let us first consider 
kfj_Il^°{k). One immediately finds from Eq. (|l3| ) that 



d'p dfjp) 
(27r)3 dp^ 



The indices l,m,n = 1,2,3 refer to the coordinates of three- vectors. The energy density carried by partons is expected 
to be finite. Therefore, /(p — oo) must vanish. Consequently, the above integral vanishes as well. Performing partial 
integration and demanding that /(p = oo) = one also proves that fc^n^™(fc) = 0. Analogously it can be also shown 
that W'^ik) = W{k). 

III. DIAGRAMMATIC APPROACH 

In this section we consider the diagrammatic Hard Loop approach to anisotropic systems. Specifically, we compute 
the QCD polarization tensor and the quark self energy for an arbitrary parton momentum distribution. 

A. Polarization tensor 

The contribution from the quark loop to the gluon self energy in the case of one quark flavor is of the form 

Kbik) -^^-^9' 1^ TTirSiqh'-Sip)] , (14) 

where q = p — k and S is the bare quark propagator. Since we are dealing with a non-equilibrium situation we adopt 
the real time formalism. Within the Keldysh representation which has been shown to be especially convenient 
in the Hard Loop approximation | [T5[ |, there are retarded (i?), advanced (A), and symmetric (F) propagators which 
in the case of massless quarks are given by 

Sr,a{p) 



p2 ± isgn(po)0+ ' 

Sf{p) = ~2m^([l ~ n(p)]e(po) + [1 - fl(p)]e(-po)) <5b') , (15) 

where n(p) (fi(p)) is, as previously, the (anti-)quark distribution function [p^ that reduces in equilibrium to the form 
(^. Performing the trace in (|lj) and suppressing the color indices, we find the retarded gluon self energy as 

2 /■ 



W'{k)^ig' J j^iq^p'' +p'^q-' ^g'^''iq-p)][AFiq)An{p) + AAiq)Ap{p)], (16) 

where Sr^a,f{p) — i> A]^^a,f{p)- Terms containing A^(<7)A^(p) and Afl(g)Afl(p) have been neglected as they vanish 
after integrating over pq. 

First, we will consider the spatial components of the polarization tensor. The other components follow from it using 
the transversality of the polarization tensor, as we will discuss below. After performing the integration over po, we 
obtain 



H'"(fc) 



g2 [ (fp fg{p) 
(2^)3 IpI 



2^'^" - k^p"" -p^k"" + (5'™(-cj|p| + k • p) 
-2w|p| + 2k-p + fc2 -isgn(|p| -uj)Q+ 

2p'p™ - fc'p™ - p'/c" + (5'™(w|p| -t- k ■ p) 
' 2w|p| + 2k • p + P - isgn(-|p| - w)0+ 



(17) 



where u) is the 0-th component of k i.e. k = (w, k) and fq{p) = n{p) + n{p) . Here the vacuum part has been neglected 
because it is suppressed compared to the matter part in the Hard Loop approximation. 

Adopting the Hard Loop approximation we assume that the internal momenta are much larger than the external 
one, i.e. w, ki <^pi. Note that for arbitrary anisotropic distributions we have to require the Hard Loop condition for 
each component of the momenta, whereas in the isotropic case w, |k| <^ |p| suffices. Expanding the expression in the 
square brackets for small external momenta yields 
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p p 



-w|p| + k • p - iO+ 



[w|p| 



k • p + ?:o+ 
k-p) 



+ 



-fc'p™ - p'fc" + (5'™(-w|p| + k • p) 



2(-tj|p|+k-p-i0+) 



2(t^|p| + k-p + iO+) 



2(-tj|p| + k • p - i0+)2 2(tj|p| + k • p + iQ+Y ■ 



(18) 



In equilibrium the first two terms vanish after integrating over p. This also holds out of equilibrium if we assume 
fqi^P) = fq{p)- Then, we arrive at the final result 



n""(fc) = -^ / d^P /(P) (fc'p'"+p'fc™)(a;|p|-k-p)+(5'™(u;lp|-k.p)2-pV'"('^'-|kn 



(2^)3 



(w|p| -k-p + i0+)2 



(19) 



where we replaced fq{p) by /(p) = n{p) + n(p) + 2Ncng{p). The point is that in the Hard Loop limit the gluonic 
contributions to the polarization tensor have the same structure as the quark ones and only the distribution 
function and the color factor chan ge. For essentially the same reason, the QCD polarization tensor, computed 
even with the complete expression ( [l8|) without assuming /(— p) = /(p), is gauge independent in the Hard Loop 
approximation. Indeed, the gluon polarization tensor has the same structure as the photon one in the Hard Loop 
limit. Since the one- loop photon polarization tensor contains no gauge boson propagator it is gauge independent. 
Consequently, the same holds for the gluon polarization tensor. The result ( pj| ) is fully equivalent to Eq. (|l^ obtained 
within the semiclassical kinetic theory. In order to show the equivalence, one performs a partial integration in (^3|) and 
immediately gets Eq. (^9|). However, we do not need to assume the reflection symmetry of the distribution function 
to derive Eq. 

Two more comments are in order here. First, if we do not assume the reflection symmetry of /(p) the first two 
terms indlq) will contribute, leading to contributions in the polarisation tensor that dominate over the Hard Loop 
result ( [l9|) and are absent in the semiclassical approximation. As explained above, the extra contribution appears to 
be gauge independent. Second, for equilibrium distribution functions (0) the integrals over |p| and over the angle in 
( |l9| ) factorize. Then, it is easy to show that (19) reduces to the well known Hard Thermal Loop result where the 
polarization tensor has only two independent components and depends on lo and |k|. 

Owing to transversality, the time-like components of H'^" follow from H'™. Indeed, H"™(fc) = A;'H''"(fc)/a; and 
H°*'(fc) = fc'fc™H'™(fc)/a;^. In order to prove the transversality of the hard loop polarization tensor in the case of 
anisotropic distributions, we compute fc^H^''(fc). Considering first the quark loop contribution, we get 



fc^H''''(fc) = z.g2 



After integrating over pq we find 



d^p 
(2^ 



[2(fc . p)p- - /A:- _ fcV] [AFiq)An{p) + A^(g)A;^(p)] 



(20) 



<r_ 

' 2 



(27r)3 



2(c^|p| 



PJIPI 



-2cj|p| + 2k • p + P - isgn(|p| - uj)0+ 
-2(-c.|p|-k.p)|p|+A;2|p| 



2aj|p| + 2k • p + A:^ — i sgn( 



CJ)0^ 



(21) 



and 



fc^H^"(fc) = 



d'p fgip) 



(27r)3 



2(u;|p| - k • p)j3" - fcV 



2^m 



-2cj|p| + 2k • p + A;2 - isgn(|p| - uj)0+ 

2(-w|p| - k ■ p)p™ - fcV" 
'2w|p| -|-2k - p + fc2 - isgn(-|p| - Lu)0+ 



(22) 



Expanding the integrands in the these expressions for small external momenta analogously to (|T^) , it is easy to show 
that ( |2l| ) and ( |2^ ) vanish in the Hard Loop approximation. This also holds if the gluon loop contribution is added as 
they have the same structure in the Hard Loop approximation. 



B. Fermion self energy 

As mentioned in the Introduction, the diagrammatic technique has the advantage over the semiclassical transport 
theory approach that it can be easily extended to fermionic self energies. Therefore, we discuss the Hard Loop quark 
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self energy for anisotropic momentum distributions. Using the Feynman gauge, the one-loop quark self energy is 
found as 



E,, (fc) = 2i Cf 5., J ^ Sip) A(q) , 



(23) 



where Cp = (N^ — 1)/Nc and now q = k — p. Adopting again the Keldysh representation, the giuon propagators in 
the Feynman gauge are given by 



g2 ± isgn(go)0+ 
Apiq) = -2ni [l + ng{q)] 5{q^) . 

Suppressing the color indces, we find for the retarded quark self energy as 



(2^ 



^ AR{p)AF{q) + Afl(p)A^(g) + AF(p)AR{q) + AAip)ARiq) 



The matter part of can be decomposed in two contributions which read after integrating over po 



(uj - p)7o - (k - p) • 7 



~2cj|p| + 2'k-p + k^ + isgTi{uj - \p\)0- 

+ Pho - (k - p) ■ 7 

2tj|p| + 2k • p + A:2 + isgn{uj + |p|)0+ 



(24) 



(25) 



(26) 



and 



(fp n(p) + n(p) 



(27r)3 



-Plo + P • 7 



-2w|p| + 2k • p + fc2 + is.g\i{uj - |p|)0- 

P7o + P • 7 

2tj|p| + 2k ■ p + fc2 + isgn(u; + |p|)0+ 



(27) 



In contrast to the polarization tensor, we need to expand the square brackets in (|2^) only to the first order for small 
external momenta, leading to 



-Plo + P • 7 



P7o + P • 7 



-2a;|p| + 2k • p - iO+ 2cj|p| + 2k • p + zO+ 



(28) 



Assuming again the reflection symmetry for the distribution functions, we obtain the final gauge independent result 
in the Hard Loop approximation 



S(fe) = ^Cf I 2rig(p) + n(p) + n(p) 
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70 + V • 7 



(27r)3 



iO+ 



(29) 



In the case of isotropic distributions (^9|) reduces to the well known Hard Thermal Loop result [g2|, where the self 
energy for massless quarks contains only two independent scalar functions depending on lo and |k|. Giving up the 
reflection symmetry of the distribution functions does not introduces new dominant terms in this case since the self 
energy follows already from the lowest order terms (p8|). We have adopted this symmetry to treat H and E in exactly 
the same way. 



IV. DISPERSION RELATIONS 



The gluon polarization tensor and quark self energy can be used to determine the dispersion relations of gluons 
(plasmons) and quarks in the quasistatic and quasihomogeneous but anisotropic state of the quark-gluons plasma. 
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A. Gluon dispersion equation 



The background gluon field ^^(fc) satisfies the following equation of motion 

k^g^'" - k^'k'' - W{k)\ A^{k) = . 



Therefore, the general plasmon dispersion equation is of the form 



det 







(30) 



Equivalently, the dispersion relations are given by the positions of the pole of the effective gluon propagator. Due to 
the transversality of H^" not all components of H^" are independent from each other and consequently the dispersion 
equation ( pO| ) can be simplified. For this purpose we introduce the color permittivity tensor e'™(fc). Because of the 
relation 

e^'^{k)E\k)E'^{k) = W{k)A^,{k)A^{k) , 
where E is the chromoelectric vector, the permittivity can be expressed through the polarization tensor as 

1 



^n'™(fc) 



(31) 



There are two other equalities which follow from the transversality of li^^^ . Namely, 

nO"(fc) = (e'™(A:) - (5'™) fc'fc" , n'°(fc) = (e'™(fc) - (5'") w ¥ 

Using the permittivity tensor the dispersion equation gets the form 



det 



2 \/m jj j^m 



with 



e'"(fc) = <5'™ + 



LU-'e'^^k) 



dfip) 







k^-v'' 



2uj J (27r)3 uj-'k-v + iO+ dp" 
In the isotropic state there are only two independent components of the permittivity tensor 

e'"(fc) = eT{k)(^5'"' ~ k^k'"/k^)^ + eL{k) /c'fc^/k^ , 

and the dispersion equation (B3) splits into two equations 



(32) 



(33) 



erik) = kVo 



eL(fc) = 



The permittivity tensor ( |33| ) was calculated for the strongly elongated parton momentum distribution /(p) and it 
was found |^,|| that there are unstable solutions of the dispersion equation (^). 



B. Quark dispersion equation 

The quark dispersion relations are determined by the poles of the hard loop resummed quark propagator or equiv- 
alently are found as solutions of the equation 



det 



^ - S(p) 







(34) 



One sees in Eq. ( p9| ) that the spinor structure of E is very simple: S(p) = 7'^S]^(p). However, we also include here 
the scalar part which is relevant for the massive quarks. Then, 



^{p)^r^^{p)+C{p) 



(35) 
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Substituting the expression (^5|) into Eq. (]3J) and computing the determinant as explained in Appendix 1 of | 
get 

[{pf^ - W{p)) {p, - E,{p)) - C^p)] ' = . 

When the momentum distribution is isotropic, the structure of S further simphfies p^ : 

Sb) = A{p) + Bip) p • 7 + C{p) . 

Then, the dispersion equation reads 

(1 - A{p)ypl - (1 - S(p))'p2 ~C^p)^0. 



, we 



(36) 



V. DISCUSSION 



In the present work we have considered an anisotropic relativistic plasma which is either in equilibrium and the 
anisotropy is caused by external fields or the plasma is out of equilibrium. In the first case we deal with the ho- 
mogeneous and static systems while in the second one it can be treated as quasihomogenous and quasistatic for 
sufficiently short space-time intervals. An example of the first case is the magnetized plasma while of the second one 
the parton system from the early stage of relativistic heavy-ion collisions where we encounter a strong anisotropy in 
the momentum distribution. 

The QCD polarization tensor has been computed in two ways. We have first applied the linear response method 
within the semiclassical transport theory and then the diagrammatic Hard Loop approach. The two methods are 
equivalent (but the distribution functions have to possess a reflection symmetry, i.e. /(— p) = /(p))- When using 
the diagrammatic approach we have referred to the real time formalism since the systems under consideration are, 
in general, out of equilibrium. According to the Hard Loop approximation, we have used bare propagators for the 
internal lines of the polarization tensor which exhibit an explicit anisotropic momentum distribution. Another method 
has been used in to study anisotropic relativistic QED plasmas in a strong magnetic field. The system has 

been assumed there to be in equilibrium and dressed propagators corresponding to electrons in Landau levels have 
been adopted. In this way, anisotropic distributions arise although the distribution functions depend only on the 
energy. 

As already mentioned, the semiclassical kinetic theory of quarks and gluons has been shown so far to be fully 
consistent with the QCD dynamics only for quasiequilibrium The considerations presented here demonstrate 

that the equivalence holds for the systems which are far from equilibrium although the space-time homogeneity must 
be invoked. Thus, the reliability of the kinetic theory methods is improved and in particular, the stability analysis of 
the parton system form the early stage of ultrarelativistic heavy- ion collisions, which has been based on the linearized 
kinetic equations, is substantiated. 

The main advantage of the diagrammatic approach over the transport one is that it allows for a systematic pertur- 
bative extension to higher order effects. Also the fermion self energy for anisotropic systems can be calculated in this 
way. Having the QCD polarization tensor and quark self energy derived here, one can construct effective gluon and 
quark propagators from the Dyson-Schwinger equation. 

The poles of the effective propagators determine (via Eqs. (|3^, ^)) the parton dispersion relations in an anisotropic 
quark-gluon plasma. In the isotropic plasma, the dispersion relations for gluons and quarks show two branches and 
start from the same energy at zero momentum | |2l| . The point is that in this case there is no direction preferred and 
the longitudinal and transverse components of the dielectric function are identical when the momentum vanishes. For 
the anisotropic systems with a preferred direction even at zero momentum, we expect additional branches and the 
degeneracy at zero momentum to be removed. In equilibrium all modes are stable or damped due to the Landau 
mechanism. In the case of anisotropic systems growing modes, i.e. instabilities, are possible. The unstable modes 
were argued to occur in the parton system from the early stage of ultrarelativistic heavy-ion collisions |^,|| . Since the 
characteristic time of instability development was estimated to be rather small (below 1 fm/c) these instabilities can 
significantly influence the temporal evolution of the parton system. 

Also the quark dispersion relations following from the effective quark propagator are of physical relevance, as they 
lead in equilibrium to interesting structures, e.g. van Hove peaks, in the dilepton production rate |p5|| , which might 
serve as a signature for the quark-gluon plasma formation. It has to be seen whether these structures also survive in 
the nonequilibrium case. 

When the plasma is in the (isotropic) equilibrium state the zero frequency limit of the longitudinal component of 
the polarization tensor in the Hard Thermal Loop limit (H^ = Hoo), which is identified with the lowest order Debye 
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screening mass, is finite. The transverse component (IIt = {Sim — kikm/\^\'^)^im/'2), on the other hand, shows no 
static magnetic screening. The situation is much more comphcated in the anisotropic plasma. The screening length 
depends on the orientation of the vector k [|j , see also |26| . 

The diagrammatic approach, following the Hard Loop resummation technique |]l0| , |l5| , allows for a systematic 
calculation of observables, such as the energy loss of energetic partons or the production of photons and dileptons 
^ . Now, the program can be extended to the anisotropic quark-gluon plasma although one has to choose a specific 
form of the parton momentum distribution. We leave all these issues for future investigations. Here we have intended 
to provide only the general formalism to study anisotropic systems of quantum fields. 
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